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Abstract. We study the Frobenius problem: given relatively prime positive integers ai, . . . , Od, find the largest value 
of t (the Frobenius number) such that X]fc=i i^kOk = t has no solution in nonnegative integers mi, . . . , ma- Based on 
empirical data, we conjecture that except for some special cases the Frobenius number can be bounded from above 

u I 5/4 

by V"i"2a3 —01—02—03. 



1 Introduction 

Given positive integers ai,...,arf with gcd(ai, . . . , o^) = 1, we call an integer t representable if 
there exist nonnegative integers mi , . . . , such that 

d 

i=i 

In this paper, we discuss the linear Diophantine problem of Frobenius: namely, find the largest 
integer which is not representable. We call this largest integer the Frobenius number g{ai, . . . , a^); 
its study was initiated in the 19th century. For d = 2, it is well known (most probably at least 
since Sylvester ^Hj) that 

g{ai,a2) = aia2 - ai - a2 (1) 

For d > 2, all attempts for explicit formulas have proved elusive. Two excellent survey papers on 
the Frobenius problem are 1 and jl5j . 

Our goal is to establish bounds for g{ai, . . . ,ad). The literature on such bounds is vast — see, for 
example, |S1 El El 13 El El • We focus on the first non-trivial case d = 3; any bound for this case 
yields a general bound, as one can easily see that g{ai, . . . ,ad) < g{ai, 02,03). All upper bounds 
in the literature are proportional to the product of two of the a^s. On the other hand, Davison 
proved in [H| the lower bound 5(01,02,03) > -v/3oi02a3 — ai — 02 — 03. Experimental data (see 
Figure below) shows that this bound is sharp in the sense that it is very often very close to 
5(01,02,03). This motivates the question whether one can establish an upper bound proportional 
to ^010203^ where p < 4/3. {p = 4/3 would be comparable to the known bounds.) In this paper 
we illustrate empirically, on the basis of more than ten thousand randomly chosen points, that 
5(01,02,03) < ^010203^/"^ - oi - 02 - 03. 

^Appeared in Experimental Mathematics 12, no. 3 (2003), 263-269. 
Keywords: The linear Diophantine problem of Frobenius, upper bounds, algorithms. 
Mathematical Reviews Subject Numbers: 05A15, 11P21; 11Y16. 
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2 Some geometric-combinatorial ingredients 



Another motivation for the search for an upper bound proportional to -^010203^ comes from the 
following formula of (21, which is the basis for our study: Let a,b,c be pairwise relatively prime 
positive integers, and define 

A''((a, 6, c) := # {(mi, 7712, ms) € : irik > 0, ami + bm2 + cm^ = t'j . 

Then 

f t / 1 1 1\ 1 /3 33 a b c 

Nt{a,b,c) = 7ri- + 7T -r + - + ]- +Ti^ - + T + - + ]- + - + -r 
Zabc 2 \ab ac be J 12 \ a b c be ac ab 



where 



+a-t{b, c; a) + cr_t(a, c; b) + a-t{a, b; c) , (2) 
{a,b;c):=- ^ 



A* 

at' ' ' 



c (A" - 1) (A'' - 1) 

is a Fourier- Dedekind sum. One interpretation of Nt{a,b,e) is the number of partitions of t with 
parts in the set {a, b, c}. Geometrically, Nt{a, b, c) enumerates integer points on the triangle 

{(xi,X2,X3) G : Xk > 0, axi + bx2 + CX3 = 1} , 

dilated by t. The Frobenius problem hence asks for the largest integer dilate of this triangle that 
contains no integer point, in other words, the largest t for which Nt{a,b,e) = 0. It is also worth 
mentioning that the condition that a, b, and c are pairwise relatively prime is no restriction, due 
to Johnson's formula if m = gcd(a, b) then 

g{a,b,e) = m g (—,—,e] + {m-l)c . (3) 

In formulas analogous to Q for d > 3 are given. In our case {d = 3), a straightforward 
calculation shows 

1 ^ g^{-2t+a+fe) 

at(a, b;e) = ; rr . 

dr ^-^ a^■n TLM^ a^■n H-MR 
k=l ~ ~ 

In fact, (Tt(a, 6; c) is a Dedekind- Rademaeher sum J^, as shown in _3_. Hence we can rewrite ((SJ as 

, , , t"^ t / 1 1 1 \ 1 /3 3 3 a b c 
Ntia, b,c) = — - + -( — + — + --)+--(- + - + - + - + — + — 
2abc 2 \ab ac be J 12 \ a c be ae ab 



, o.n ^ sin ^ ~ 46 sin ^ sin ^ ~ ^ 

k=l a a k=l b k=l 



4:a sin — sin — 46 ^ sin ^ sin ^ 4c ; sm — sm 
If we write the "periodic part" of Nt{a, b, e) as 



1 ^ g2lfc(2t+fe+c) ^ ^ ^^(2t+a+c) ^ ^Zlfc(2t+a+fe) 

Pt{a, 6, c) .- — 2^ sin ^ sin ^ ^ 46 ^ sin 4^ sin 4^ ^ 4^ sin ^ sin ^ 

fc=l a a k=l o k=l c c 
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@ becomes 

If we can bound Pt{a, b, c) from above by, say, B then the roots of Nt{a, b, c) — and hence g{a, b, c) — 
can be bounded from above: 

b, c) < a6c (-i + + ^) + ^ \ He + + - ih ih + I + I + Tc + + f,) - B)) 

= - U<^ + b + c) + ^Uabc)m ; + ^^ + l)'-labc{l + l + l + ^^ + l^ + ^)+2Babc 

= ^2B a6c+^(a2 + 62 + c2) -\{a + b + c) 

From this computation, the question of the existence of an upper bound for g{a, b, c) proportional to 
V abc comes up naturahy. Unfortunately, it is not clear how to bound the periodic part Pt{a, b, c) 
effectively. An almost trivial bound for Pt{a,b,c) yielded in |3j the inequality 

g{a, b,c) < - ^a/ abc {a + b + c) — a — b — cj , 

which is of comparable size to the other upper bounds for g{a, b, c) in the literature. However, we 
believe one can obtain bounds of smaller magnitude. 



3 Special cases 



On the path to such "better" bounds, we first have to exclude some cases which definitely yield 
Frobenius numbers of size a^. One of these cases are triples (a, b, c) such that c is representable by 
a and b: by we obtain in this case g{a, b,c) = ab — a — b. 

A second case of triples (o, b, c) that we need to exclude are those for which a\{b + c). Brauer and 
Shockley ;5. proved that in this case 



g{a, b, c) = max I b 



ac 
b + c 



ab 
b + c 



a . 



Here [xj denotes the greatest integer not exceeding x. 

An even less trivial example of special cases was given by Lewin jTT], who studied the Frobenius 
number of almost arithmetic sequences: If m,n> 0, gcd(a, n) = 1, and d < a, then 



g (a, ma + n, ma + 2n, . . . , ma + (d — 1) n) 



m 



d-1 



+ m — l\a+{a — l)n 



For arithmetic sequences (m = 1), this formula goes back to Roberts consecutive numbers 

(m = n = 1) it is due to Brauer j^. For the special case d = 3, we obtain 



g{a, ma + n, ma + 2n) 



m 



a 
12 



1 ) a + (a — 1) n . 
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As a function in a, b := ma + n, c := ma + 2n, this Frobenius number grows proportionally to 
ab, which means an upper bound proportional to \fabc with p < 4/3 can not be achieved. Hence 
in our computations and conjectures about upper bounds for g{a,b,c), we will exclude the cases 
of one of the numbers being representable by the other two, one number dividing the sum of the 
other two, and almost arithmetic sequences. Finally, as noted above, thanks to © we may assume 
without loss of generality that a, 6, and c are pairwise coprime. The triples (a, b, c) that are not 
excluded will be called admissible. 



4 Computations 

In the present section we discuss the computation of the Frobenius number. For convenience we 
computed the number /(a, b, c) = g{a, b, c) + a + b + c. It is not hard to see that /(a, b, c) is the 
largest integer that can not be represented by a linear combination of a, b, and c with positive 
integer coefficients. The respective counting function 

A^t(a, ft, c) := # {(mi, m2, ms) G : mk > 0, ami + bm2 + cms = t} ■ 

can also be found in 3 and is closely related to Nt{a, b, c): 

— ^ / I I i\ 1 /3 33 a b c 

Nt{a,b,c) = -— --(- + - + -) + —(- + - + - + - + - + - 
Zabc I \ab ac be J 12 \a b c be ae ab 



Trkc AU / J ■ Trka 



4a ^ ^in ^ 'iin ^ 4b ^ tjin ^ 'jin ^ 4c ^ 'jin ^ sjin ^ 

sm ^ sm ^ ^0 sm ^ sm ^ sm ^ sm ^ 

The following illustrates our algorithm. 

STEP 0: Initiate the intervals II, 12, 13 for the selection of the 
arguments a,b,c; 

STEP 1: Draw at random integers a,b,c from II, 12, 13 respectively; 

STEP 2: Test a,b,c for coprimality and for almost arithmetic sequences; 

STEP 3: IF (a,b,c are not pairwise coprime) or IF (a,b,c are almost arithmetic) 
{discard a,b,c and GOTO STEP 1} 
ELSE {SET delta <- min(a,b,c); GOTO STEP 4>; 

STEP 4: Compute z=sqrt (a*b*c) , 

SET mb <- INT(sqrt(3)*z)+delta; 
SET t <- mb; 



STEP 5: Compute NB(t,a,b,c); 
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STEP 6: IF (NB (t , a,b , c) >0) {SET t <- t-1, and GOTO STEP 5> 
ELSE {GOTO STEP 7>; 

STEP 7: SET f <- t; 

STEP 8: IF(mb-f < delta) {SET 
mb <- mb+delta 
t <- mb 
GOTO STEP 5> 
ELSE {GOTO STEP 9>; 

STEP 9: PRINT f(a,b,c) <- f; 
STOP. 

For example, for a = 7,b = 13, c = 30 the program yields the Frobenius number /(7, 13, 30) = 95, 
or g{7, 13, 30) = 45. This program was tested against arguments which yield known results, and 
found to be correct. 

Our program is to choose at random arguments a, b,c in a certain range (in our case [1, 750]), and 
test the triplets for admissibility. For admissible triplets a, b, c we compute the Frobenius number 
f{a,b,c) based on the straightforward observation that, once we have a = mm(a,b,c) consecutive 
integers which are representable, we know that every integer beyond that interval is representable 
as well. We start searching for roots of Nt{a, b, c) at the lower bound \/3abc. If a root is found at an 
integer /, we repeat this search until we found an interval of a integers t with Nt{a, b, c) > 0, that 
is, an interval of a representable integers. At this stopping point, the integer / is the sought-after 
Frobenius number f{a,b,c). 

We have created a PARI-GP program^, following the above algorithm. The program proved to be 
quite efficient, since most of the values of f{a,b,c) were found to be close to the lower bound 
V^abc, as shown in the analysis below. The Dedekind-Rademacher sums appearing in (jlj) can be 
computed very efficiently because they satisfy a reciprocity law (^Sji for computational complexity 
see also ITUI), which allows us to calculate their values similar in spirit to the Euclidean algorithm. 
This implies that for a given t, Nt{a,b,c) can be computed with our rather simple algorithm in 
0(log(c)) time, assuming that c = max(a, 6, c). Hence if f{a,b,c) is close to the lower bound 
VSabc — which, again, happens in the vast majority of cases — , we obtain f{a,b,c) in 0(alog(c)) 
time. On the other hand, we can of course not assume that f(a,b,c) is close to VSabc; still we 
get 'at worst' a computation time of 0(a61og(c)). What makes this analysis even more appealing 
is that it applies to the general case of the Frobenius problem. As mentioned above, there is an 
analog for ^ and @ for d > 3 [3], which again is a lattice-point count in a polytope and as such 
is known (for fixed d) to be computable in O {p (logai, . . . , loga^)) time for some polynomial p 2 . 
With an analogous algorithm for the general case, we would hence be able to compute /(ai, . . . , a^) 
in O (aia2P (logoi, . . . , log Orf)) time, where oi < 02 < • • • < a^. As in the three- variable case — 
in fact, even more so — , most Frobenius numbers will be situated very close to the lower bound 

^Our program can be downloaded at www. math. binghamton.edu/matthias/f robcomp.html 
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■v/SoioifflS) which means that in the vast majority of cases we can expect a computation time of 
O (aip(logai, . . . ,log ad)). 

The computational complexity of the Frobenius problem is very interesting and still gives rise to 
ongoing studies. Davison [HI provided an algorithm for the three-variable case {a < b < c) which 
runs in 0(log5) time. The general case is still open. While Kannan [^j proved that there is a 
polynomial-time algorithm (polynomial in log ai, . . . , log a^) to find g{ai, . . . ,ad) for fixed d, no 
such algorithm is known for d > 3. The fastest general algorithm which we are aware of is due to 
Nijenhuis jl2j and runs in 0{daloga) time, where a = min(ai, . . . ,0^). Hence, while our primitive 
algorithm is not competitive for the three-variable case of the Frobenius problem, it might be 
worthwhile to develop it further in the general case. 

We initially implemented our program as an MS-DOS QUICK BASIC program and experienced some 
interesting problems due to floating-point errors: computing generalized Dedekind sums can get 
challenging for large arguments. These problems were only discovered when we reimplemented the 
algorithm in PARI-GP, which has an extended precision aritmetic and also keeps track of roundoff 
errors effectively. It is worth mentioning that both Knuth's algorithm [10, for the computation 
of Rademacher-Dedekind sums and Davison's algorithm 6 for computing g{a, b, c) are integer 
algorithms and therefore are very stable. 

With our program we generated at random 10000 admissible triplets. Our main question is the re- 
lation of the Frobenius number / = /(a, b, c) to z := V abc. The following is a statistical description 
of the ratios R := f /z. 

Descriptive Statistics 

Variable N Mean Median StDev Min Max Qi Q3 
R 10000 2.283 2.012 0.737 1.736 9.332 1.940 2.299 

Qi and Q3 are the first and third quartiles, respectively. We see in the above table that 50% of the 
cases have a ratio smaller than 2.01, and 75% have ratio smaller than 2.30. In the following figures 
we present a box-plot and a histogram of the variable R. 

In the box-plot, the bottom line of the box corresponds to the first quartile Qi. The top line of the 
box corresponds to the third quartile Q3. There are 980 points above the value of i? = 3. Only 24 
points, which are listed below, have a value of R greater than 7.5. 
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Figure 1: Box-plot 
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In Figure 01 we presents all the points {z, /). Notice that all the points are above the straight line 
^/Sz, which illustrates Davison's lower bound. The upper bound is, however, convex. It is included 
in the figure as the graph z^^^. 



5 Conjectures and closing remarks 



Randomly chosen admissible arguments tend to yield a Frobenius number / smaller than the 
expected number (mean) which is estimated to be 2.282. The distribution of i? = f/z is very 
skewed (positive asymmetry) as seen in the Figure ^ Since 10000 random points yielded / < z^/^, 

or g{a, b,c) < v abc — a — b — c, the probability that a future randomly chosen admissible triplet 
with z = Vabc < 20000 will yield / > z^/^ is smaller than 1/10000. 

In general, our data suggests that one can obtain an upper bound of smaller magnitude than what 
the above cited results state. Again, the upper bounds in the literature are comparable to an upper 

4/3 

bound proportional to v abc . We believe the following is true. 

Conjecture 1 There exists an upper bound for {a,b,c) proportional to V abc' where p < |, valid 
for all admissible triplets {a,b,c). 

In fact, our data suggests, more precisely, that for all admissible triplets (a, 6, c), 

/ 5/4 

g{a, b,c) < V abc — a — b — c . 
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Figure 3: / = /(a, b, c) as a function oi z = vabc 



It is very improbable that a randomly chosen admissible triple (o, 6, c), such that Vabc < 20000, 

will yield g{a, b,c) > V abc —a — b—c. However, we remark that there might be specific structures 

of triples {a,b,c), close to almost arithmetic, for which g[a,b,c) > y abc —a — b — c. This is 
generally not the case. 
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